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SYSTEM FOR DUST DETECTION AND DUST IMAGE 
ENHANCEMENT 

[00011 BACKGROUND OF THE INVENTION 
[0002] Field of the Invention 

[0003] This invention deals generally with digital enhancement of meteorological 
imagery of dust clouds over ocean backgrounds. 

[0004] Description of the Related Prior Art 

[0005] The detection and tracking of airborne dust from satellite has been a long-standing 
interest of the Navy and academia alike, and there is no shortage of papers treating the 
subject (e.g., Shenk, W. E.. and R. J. Curran. 1974: The detection of dust storms over land and 
water with satellite visible and infrared measurements, Mon. Weather Rev., 102. 830-837. 
Carlson. T. N.. 1978: Atmospheric turbidity in Saharan dust outbreaks as determined by analyses 
of satellite brightness data. Mon. Weather Rev.. 107. 322-335. Legrand. M.. M. Desbois. and K. 
Vovor. 1987: Satellite detection of Saharan dust: optimized imaging during nighttime, J. 
Climate. 1. 256-264. Jankowiak. I., and D. Tanre. 1992: Satellite climatology of Saharan dust 
outbreaks: method and preliminary results, J. Climate. 5, 646-656.) 

[0006] These earlier works made use of contrast and spectral signatures to identify and 
retrieve properties of dust against land and ocean backgrounds. 

[0007] A recent study (R. Cantu; "The Role of Weather in Major Naval Aviation 
Mishaps", MS Thesis, 106 pp., NPS Monterey) has indicated that poor visibility has 
factored into approximately 54% of Navy aircraft "Class-A" mishaps, for an estimated 
annual loss of S51 million. Of these mishaps it was reported that 56% were preventable 
had better forecasting and/or observational tools been available. The current Dust 
Enhancement Product represents a paradigm shift in the way dust is observed from the 
space platform, and is directly applicable to mitigating this multi-million dollar liability. 
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[0008] Alternative methods exist for enhancing dust over land and water. Single-channel 
imagery with scaling thresholds chosen to enhance a small dynamic range of temperature 
or reflectance values has been applied to visible and infrared geostationary imagery. The 
resultant enhancement, usually only of marginal visual quality, is highly dependent on the 
thresholds chosen and hence not suitable for operations where optimal values for these 
thresholds may change dynamically over space and time. Combinations of single- 
channel visible and infrared data is a marked improvement over the scaling described 
above, but still suffers cloud/dust ambiguity for lack of multi-channel visible data (i.e., 
the basis for the over-water Dust Enhancement Product, which requires blue and short 
wave infrared channels). The Empirical Orthogonal Function (EOF. or sometimes 
referred to as " principle component") method has been shown to do a good job of 
enhancing dust over land and water. However, results produced by this method again are 
highly dependent on the variability of the scene and hence will provide an inconsistent 
enhancement depending on background and dimension of the region. Spatial contrast 
reduction, which takes advantage of the blurring affect of dust over otherwise sharply- 
defined terrestrial features, has been applied successfully for dust detection over land, but 
naturally fails over ocean and is limited in regions where terrain is less variable (e.g., 
desert plains, where dust incidentally is most common). 

|0009] None of the above-mentioned methods enlists multi-channel visible data to take 
advantage of the inherent spectral variability of dust at these wavelengths. The unified 
land/ocean Dust Enhancement Product is a unique solution to the daytime dust detection 
problem. A comparison between the visible/infrared combination technique (applied 
here to Meteosat 5 data) and the current algorithm (applied to Terra-MODIS data) for a 
space and time co-located dust storm event in Southwest Asia has been conducted. The 
superior detection of dust by the current method over both land and water as 
demonstrated by these data are thought to be of immediate relevance to Navy 
Meteorology/Oceanography (METOC) operations in several regions of the world that 
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experience significant dust outbreaks (including Southwest Asia, Africa, the 
Mediterranean, and the West Pacific). 

[0010] SUMMARY OF THE INVENTION 

(0011] The dust enhancement system of the invention is to provide an indication of the 
horizontal distribution of airborne dust over both land and water surfaces for daytime 
scenes. The regions of dust are enhanced as shades of pink, depending on opacity, 
against darker non-dust backgrounds. 

[001 2] The primary advantage of the Dust Enhancement Product is its ability to call 
immediate attention to areas where there airborne dust resides, often in regions where 
even the true color imagery experiences difficulty in discerning. The enhancement 
highlights mesoscale circulations in optically thin dust, and is also useful for identifying 
volcanic plumes (e.g., Mt. Etna 7/26/01 eruption— image attached). The latter capability 
is of particular relevance to both military and civilian aviation. The algorithm is 
applicable to a wide range of ocean color instruments in the ocean-only mode (i.e., true 
color over land surfaces) or in the ocean/land mode provided infrared information is 
available either on the same platform or co-registered with temporally matched imagery 
from an independent observing system. 

|0013] The techniques of the invention demonstrate that the current technique, as applied 
here to the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) radiometer, has 
significant potential to benefit the warfighter through its ability to provide through high 
spatial resolution satellite imagery a detailed depiction of atmospheric dust outbreaks 
with enhancements for optically thin regimes, as well as other meteorological/land- 
surface features of interest (e.g., clouds, fires, snow cover). Based on these findings, 
NRL and Fleet Numerical Meteorology and Oceanography Center (FNMOC) have 
collaborated to host a subset of these products upon the FNMOC Secure Internet 
(SIPRNET) beta-webpage for the purpose of i) making this new technology available to 
Navy assets, and ii) receiving feedback useful for improving upon its current utility. 
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10014] BRIEF DESCRIPTION OF THE DRAWINGS 

10015] Fig- 1 illustrates a SeaWiFS sensor (left) flying aboard the SeaStar platform 
(right). 

|0016] Fig. 2 illustrates raw SeaWiFS data from a Navy Regional Center (NRC) 
receiving station (left) and the noise/line corrected equivalent (right). 

[0017] Fig. 3 illustrates a SeaWiFS High Resolution Picture Transmission (HRPT) true 
color image of the Eastern Atlantic Ocean, Iberian Peninsula, and Northwestern Africa, 
captured by NRC in Rota, Spain (02/13/01, 1255Z). A large dust plume wrapping into a 
baroclinic system is indicated. 

[0018] Fig. 4 illustrates a SeaWiFS HRPT true color image of the Persian Gulf area, 
captured by NRC in Bahrain. Airborne dust is indicated. Phytoplankton is responsible 
for the green pigmentation in the littoral regions. 

|0019] Fig. 5 illustrates a SeaWiFS HRPT true color image of the Korean Peninsula, 
captured by NRC in Yokosuka Japan. Yellow arrows indicate areas of biomass 
smoke/pollution. 

|0020] Fig. 6 illustrates a true color (top) and vegetation enhancement (bottom). Bright 
red indicates regions with abundant plant life. Note seasonal changes in Fall vegetation 
compared against Summer (Fig. 5). 

[0021 J Fig.7 illustrates preferential absorption of blue light by mineral dust, resulting in 
yellow reflected light. Clouds scatter more uniformly over the visible wavelengths, 
resulting in white light. 

[0022] Fig. 8 illustrates examples of thin dust enhancements over the waters of the 
Persian Gulf. Optically thin dust plumes appear as darker tones of red. Land regions in 
the dust enhancement have been masked. 

[0023] Fig. 9 illustrates true color imagery of Typhoon Pabuk making landfall on Japan. 
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[0024] Fig. 10 illustrates an enhancement captured volcanic ash plume from the 2001 
eruption of Mt. Etna. Land regions in the dust enhancement have been masked out. 

[0025] Fig. 11 illustrates a MODIS dust enhancement (0.5 km resolution) over the 
Northern Arabian Sea, true color over land. Dust over the ocean and appear as shades of 
pink. Dry lakebeds over land appear as patches of pink. Accompanying text from Navy 
MetOc post deployment report corresponding to the storm observed in this image. 

[0026] DETAILED DESCRIPTION OF THE INVENTION 

[0027] The Satellite Meteorological Applications Section at the Naval Research 
Laboratory (NRL) in Monterey has developed true color, vegetation enhancement, and 
airborne sand/dust enhancement processing capabilities for the Sea-viewing Wide Field 
of View Sensor (SeaWiFS) instrument. These products are made in near real time 
(typically 60 to 90 minute latency from time of ingest at the receiving stations) using 
telemetries captured at Navy Regional Centers (NRCs) located in Spain (Rota), Bahrain, 
and Japan (Yokosuka), and populate a password-protected website for several predefined 
sectors within the general coverage regions of these centers. Presented here is an 
overview of the SeaWiFS sensor and its capabilities/limitations, a summary of the 
procedure, and the current processing architecture. Examples drawn from the archived 
database (automated processing of these data commenced August 8, 2001) are presented 
to illustrate their quality and potential operational utility. 

[0028 1 Platform/Instrument Specifications 

[0029] The SeaWiFS instrument (Firestone and Hooker, 1996), launched aboard an 
extended Pegasus rocket on August 1, 1997, is an 8-channel whiskbroom scanning 
radiometer featuring 8 narrowband (-20 nanometers (nm) wide) spectral channels 
centered at 412, 443, 490, 510. 555, 670, 765. and 865 nm wavelengths. Carried aboard 
the SeaStar spacecraft (see Fig. 1), SeaWiFS follows a sun-synchronous (98.3° 
inclination. 705 kilometer (km) altitude above mean sea level (AMSL)) orbit with a local 
noon descending equatorial crossing time. It has a local-area-coverage (LAC) swath 
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width of 2801 km (corresponding to a maximum scan angle of 58.3°) and a nominal 
spatial resolution of 1.1 km at nadir. Global Area Coverage (GAC) is also available at 4 
km degraded resolution and a reduced 1502 km swath (data limited to 45° scan angle). 
The 10-bit digitized data are transmitted down to Earth at an L-band frequency of 
1.70256 Gigahertz (GHz). The instrument scans across 1285 picture elements at a rate of 
6 Hz. Typical High Resolution Picture Transmission (HRPT) station telemetries include 
typically 4000 to 5000 scan lines of data (depending on local 
topographical/anthropogenic obstructions), translating to 11-13 minutes of capture time 
from the nominal initial-capture time stamp on the filename. 

[0030] The primary function of SeaWiFS is to detect and monitor ocean color. 
Specifically, it has been designed to retrieve the small component (-5% of total signal) of 
water leaving reflected radiance whose spectral signature provides useful information 
about ocean phytoplankton content (expressed in terms of chlorophyll-A retrievals). A 
varying along-track sensor tilt angle of +/- 19.8° is used to avoid regions of sun glint 
(specular reflection of the solar disk upon the ocean surface) where chlorophyll-A 
retrievals cannot be performed. The information obtained from such retrievals is 
valuable to marine research and industry interests. Applications include guidance to the 
fishing industry, diver visibility information, surface currents, physical oceanography, 
general marine biology, and climate impacts on the marine ecosystem. While 
atmospheric constituents such as cloud, dust, and aerosol are considered as noise to ocean 
color algorithms, the current research regards these components as the signal of interest. 
This description of the invention report details how SeaWiFS data are applied to enhance 
atmospheric and surface features/phenomena. 

[0031] Instrument Advantages 

[0032] The main advantage provided by SeaWiFS over optical spectnim instruments 
flown on the geostationary (e.g., the Geostationary Operational Environmental Satellite 
(GOES) constellation, Meteosat 5,7, Geostationary Meteorological Satellite (GMS), and 
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the Indian National Satellite System (INSAT)) and other polar orbiter (e.g.. National 
Oceanographic and Atmospheric Administration (NOAA), Defense Meteorological 
Satellite Program (DMSP), and many others) platforms for the current application is its 
multiple narrowband visible channels that provide a true color imaging capability. 

[0033] Instrument Limitations and Other Caveats 

[0034] The main deficiencies of SeaWiFS in terms of 24-hour operational utility are its 
lack of infrared channels, its limited temporal resolution, and proprietary nature. Infrared 
data are useful for nighttime imaging, cloud particle size estimation, and airborne dust 
enhancement over land surfaces, among other uses. At mid to low latitudes, only one to 
three passes are available near local noon — providing at best a "snapshot" of the 
synoptic-scale situation on any given day. Captures are increasingly frequent at the 
higher latitudes owing to the orbital geometry. Typical captures for the regional centers 
currently processed are on the order of 2 to 3 per day. Thus, SeaWiFS products are most 
useful when juxtaposed with geostationary data that offer the temporal information 
required to infer details of the atmospheric circulation. 

[0035] Because SeaWiFS was built, launched, and operated by Orbital Sciences 
Corporation (OSC), the data are formally available to the general public on a research- 
only basis and are provided at cost. The National Aeronautics and Space Administration 
(NASA) has contracted with OSC for the research use of a 5-year data set commencing in 
September 1997. Researchers who have registered as authorized SeaWiFS user are 
granted permission to obtain the data from the NASA SeaWiFS archive, approved NASA 
Space Act Agreement ground station, or authorized individual. As authorized receiving 
stations, the NRCs providing near-real-time LAC data for this research are covered under 
the existing agreement, with the only criterion being that NRL does not disseminate the 
LAC-derived products to the general public in near real-time. 

[0036] Data Processing 
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|0037] Several options exist for processing the SeaWiFS telemetry, depending on the 
format of the data. The formats currently received at NRL are NASA Level-0 (L0; 10-bit 
Data Capture Facility frame-formatter files), with Terascan Data Format (TDF) (pre- 
calibrated and navigated) upon special request. Data from the NRCs are in NASA L0 
format, and require several preprocessing steps to calibrate and navigate the data as 
described below. For reasons described further along, the current approach to creating 
SeaWiFS products is based on NASA L0 data and is independent of the Terascan 
software package. Several approaches have been iterated upon in arriving at the current 
procedure including various blended combinations of software packages — some more 
appropriate for rapid deployment to remote operational centers and others geared toward 
in-house processing with all the conveniences of locally available third party software. 
While the purpose of this report is to summarize the current approach, alternative options 
and their associated data/software requirements will be touched upon in passing. 

[0038] Worth noting is that the Perl programming language, running on the Linux 
platform, was selected as the scripting language for automated processing of SeaWiFS 
data. Perl draws from the strengths of many popular scientific programming languages 
such as C++ and Fortran while enhancing on the capabilities of other standard shell 
scripting languages. The Perl scripts written in support of the current SeaWiFS 
processing perform the requisite secretarial operations of data retrieval, file and directory 
manipulation, execution and status monitoring of telemetry updating and data processing 
codes, dissemination of imagery products to the website, and general housekeeping. Perl 
remains the lone constant among all processing options discussed in this report. 

[0039] SeaDAS Processing Option 

|0040] One of the available SeaWiFS processing options is the GSFC SeaWiFS Data 
Analysis System (SeaDAS) software package (Fu et «/., [1998]), which was designed 
specifically for viewing and basic manipulation of SeaWiFS data. It uses an embedded 
runtime Interactive Data Language (IDL) license as its operating paradigm. The SeaDAS 
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software is designed to operate on NASA L0, Level-IA (L1A; raw spacecraft and 
instrument telemetry retained with convened geo-location, instrument telemetry, and 
selected spacecraft telemetry), and Level-IB (LIB; sensor calibration applied to the 
data). LI A and LIB data are stored in Hierarchical Data Format (HDF). The software 
includes a graphical user interface enabling straightforward viewing and manipulation of 
the data, as well as command-line functions useful for incorporating into scripts as batch 
jobs. The SeaDAS package was used early on in this research to examine the general 
utility of SeaWiFS data in terms of coverage and data quality before proceeding with 
customized application. While SeaDAS is not used in the current processing (in the 
interest of gaining a higher level of control over product display options), several 
functions were adopted from this package for conversion from NASA L0 to LIB HDF 
files and data correction. SeaDAS remains an extremely useful stand-alone tool for 
SeaWiFS data visualization. Versions of SeaDAS for several operating platforms are 
available and a limited level of helpful online support exists. 

1 0041] The Currently Implemented Processing Scheme 

[0042] A simple overview of the current processing flow is as follows: I) SeaWiFS LAC 
telemetry captured at an NRC site and converted immediately to NASA L0 fonnatted 
file, II) Data sent to the satellite data storage facility at Naval Oceanographic Center 
(NAVO) at Stennis Space Center in Mississippi via file transfer protocol (ftp) push, III) 
Data are pushed via ftp to NRL Monterey from NAVO, IV) A periodic local cron job 
searches for receipt of new files, and if found: V) SeaWiFS processing script commences 
to calibrate, clean, apply atmospheric corrections, and generate products, VI) Products 
are uploaded to a password protected NIPRNET website. Latency of near real-time 
products is generally on the order of 1 .5 hours from time of ingest at the NRC site. The 
remainder of this section provides additional details on the steps outlined above. 

10043] NRC Stations 
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[0044] The mid-latitude NRCs included in the current NRL processing are Rota [36.63N. 
6.34E], Bahrain: [26.16N, 50.39E], and Yokosuka [35.28N, 139.67E]. Each station 
captures SeaWiFS passes on or around its local noontime hours (1 100 — 1400 Z for Rota, 
0700—1000 Z for Bahrain, and 0130—0430 Z for Yokosuka). NRL receives the NASA 
L0 data from NAVO as Gnu-compressed files typically 1 hour after the nominal capture 
time. File sizes vary according to pass coverage, with typical values ranging from 20 to 
80 Megabytes. The 55 Gigabyte hard drive currently used as a storage bay for raw data, 
algorithms, and products is sufficient for several months of online storage. All data are 
eventually archived to CDROM or tape storage in a hierarchical structure partitioned by 
station, year, month and day. 

|0045] Data Calibration 

|0046] A cron job launched every half hour initiates the processing of newly received 
SeaWiFS data. A loop over all files present in the ftp "receive directory" checks for any 
new and unprocessed data. If a new file has arrived, it is uncompressed and copied to the 
local storage drive for subsequent processing to LIB calibrated radiances (expressed in 
mW/cm 2 -sr-um). Binary executables taken from the SeaDAS software package are used 
for this purpose, and require an updated SeaStar ephemeris for correct earth location 
(SeaWiFS requirement of 1 -pixel geolocation accuracy). These data are available and are 
updated on the local system if the current elements file is found to be more than one day 
old. Any errors incurred during the uncompressing or calibration steps are reported to a 
log file and the script is exited gracefully. 

1 0047] Sector Definition and Noise Corrections 

[0048] Upon completion of LIB processing, an IDL script for customization of the data 
is spawned. Latitude, longitude, sun/sensor geometry, and SeaWiFS channel variables 
are extracted from the HDF data. A table of pre-defined sectors for each station is 
included and readily updated for new operational areas of interest or one-time case 
studies. Definition of these sectors is based on the specification of latitude/longitude 
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boxes (lower left and upper right corners). Owing to the scanning geometry, regions near 
the edge of the swath suffer spatial resolution degradation. As a result, imagery for these 
pre-defined sectors may include only partial coverage at compromised quality. A 
minimum data check of 20% coverage over the sector in question is enforced prior to 
proceeding with the processing. To ensure that the optimal information content for 
every SeaWiFS pass is provided, a "floater" sector is included which selects a box data 
away from the swath edges (poor resolution) and the horizons (avoiding numerous scan 
line drop-outs). The trade off for high image quality is a random coverage area from 
pass to pass. 

[0049] Once the sector of interest has been defined (or loops over numerous sectors), an 
embedded loop over SeaWiFS channels is performed for corrections of pixel noise and 
line dropouts. Noisy pixels appear in the imagery as a "salt and pepper effect" owing to 
the presence of spurious very high or very low data values. Line dropouts appear in 
imagery as black scan lines devoid of data or filled with noise. These imperfections are 
present in the data due to a combination of possible hardware limitations/problems and 
clutter associated with nearby structures. Cleaning of the data is accomplished in two 
ways. Lines identified as dropouts are replaced by previous good scan lines (bad lines at 
the very beginning or end of a data stream are replaced by their nearest good neighbor). 
If multiple adjacent dropouts are present, the previously patched line is used once again. 
This correction loses its cosmetic effectiveness for large data gaps, but provides marked 
improvements in cases of spurious line dropouts. Noisy pixels, assumed to be spurious, 
are identified using a threshold based on the absolute values of pixel derivatives 
computed along a scan line. When flagged, they are replaced by an average of the two 
adjacent pixels on the same scan line. Noise pixels on the edges of lines are simply 
replaced by the value of their nearest neighbor. Fig. 2 illustrates the result of corrections 
applied to a subsector of SeaWiFS data collected over the Alps. The overall effect of 
these bad line and pixel noise corrections is dramatic, and regarded as an essential step in 
the production of high-quality SeaWiFS imagery. Additional noise corrections (not yet 
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implemented) may be required to correct pixels where one of the channels is missing data 
(but was not flagged by the noise threshold). This results in a magenta, yellow, or cyan- 
biased pixel values in the RGB products. 

[0050] Atmospheric Corrections 

[0051] The effects of atmospheric molecular (sometimes referred to as "'Rayleigh") 
scattering must be removed from the data before creation of true color products. The 
amount of Rayleigh scatter is proportional to X" 4 (wavelength) and hence is strongest for 
the shorter (e.g., blue, violet) wavelengths. If left uncorrected it will manifest in the 
imagery as a semi-transparent milky blue haze increasing in opacity toward the edges of 
the swath (as higher sensor angles peer through an optically thicker atmospheric path), 
detracting from image contrast and clarity. To remove this undesirable component of the 
satellite signal, we use a radiative transfer model (Miller et al, [2000]) to simulate the 
clear sky Rayleigh scatter over a dark surface as a function of the solar zenith, satellite 
zenith, and solar-satellite relative azimuth angles. Pre-calculated look-up tables were 
created for each SeaWiFS channel and the correction was applied on a pixel-by-pixel 
basis, provided the sun/satellite geometry in the Level IB data. The primary caveat 
associated with this correction is the assumption of a cloud-free Rayleigh contribution 
that in some instances may result in an over-correction of pixels containing high clouds 
near the edges of the swath. Visual inspection of imagery has not indicated a significant 
distortion owing to this effect, however. Regions near the edge of the SeaWiFS swath 
should in general be interpreted with cautious regard to the possible distortions owing to 
larger atmospheric path lengths, three-dimensional effects (parallax), and possible bow- 
tie effects. 

[0052] Scaling 

[0053] Once the data have been calibrated, cleaned of noise, and corrected for Rayleigh 
scatter, the spectral radiances are converted to equivalent isotropic reflectances (Rx) 
according to the relationship: 
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(1) 



where u 0 is the cosine of the solar zenith angle and F 0 is the solar spectral flux at 
wavelength X. By definition, the isotropic reflectance represents the reflectance produced 
by an equivalent lambertian surface, and varies from 0.0 to 1.0. Any spurious values 
greater than 1.0 or less than 0.0 are thresholded by these upper and lower bounds. The 
reflectances are scaled by log\o in an attempt to equalize the imagery and prevent the very 
bright cloud/snow pixels from dominating the information content of the products. These 
scaled reflectances are byte-scaled over the 256-element color table range of the IDL 
palette, with a minimum logarithmic reflectance of -1.65 mapped to the minimum 
saturation value (0) and 0.0 mapped to the maximum (255). 

[0054] The scaled data are warped to a mercator projection corresponding to the sector of 
interest. This procedure is currently the most computationally time-consuming 
component of the processing, and options for optimization or alternative approaches are 
currently being explored. 

[0055] The final image products are created as jpegs at 80% of full quality. File sizes 
range from 200 to 800 Kilobytes depending on the sector and image information content. 
The true color, vegetation enhancement, and dust enhancement products for various 
stations and subsectors are copied to a password-protected online holding directory on the 
NRL Satellite Meteorology web page. 

[0056] True Color Enhancement 

|0057] To produce the true color images, the blue (412 nm), green (555 nm), and red 
(670 nm) SeaWiFS data are loaded into the respective blue, green, and red IDL "color 
guns." The data comprise a single matrix of dimension (n.m,3) where n=along-scan 
direction and m=along-track direction. The resultant imagery is similar, but not identical, 
to what would be observed by the human eye from the vantage point of outer space: 
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clouds as shades of gray to white, ocean surfaces as shades of blue, and land features as 
earth-tones and green in regions of heavier vegetation. 

[0058] Examples of true color SeaWiFS imagery are shown for the Rota. Bahrain, and 
Yokosuka regional centers in Figs. 3-5, respectively. Fig. 3 demonstrates the immediate 
ease with which significant airborne dust is identified over oceans. The image depicts a 
widespread easterly plume of Saharan dust captured within the synoptic flow of a 
baroclinic system off the coast of Spain and Portugal. SeaWiFS is useful in monitoring 
these dust plumes as they migrate across the Atlantic Ocean (often reaching the eastern 
United States). The green colors in littoral regions of the Persian Gulf in Fig. 4 
correspond to a higher amount of bioactivity (e.g., plankton) in nutrient-rich waters. 
Cloud streets over land and evidence of biomass burning in the form of diffuse smoke 
plumes are apparent over Southeastern Asia in Fig. 5. The NRC in Yokosuka also 
captures seasonal dust storms blowing off the Gobi desert. Owing to the general 
circulation of the atmosphere, large Gobi and Saharan dust storms may reach the 
continental United States (the former crossing the Pacific with the "Roaring Fourties", 
and the latter crossing the Atlantic with the Tropical Trades). 

[0059] Vegetation Enhancement 

(0060] Because green vegetation possesses a relatively high albedo compared to soil in 
the "reflective" or shortwave infrared (0.7-1.3 jam), a combined vis/shortwave-IR product 
that enhances vegetation regions is generated in a way similar to the true color 
enhancement described above. This vegetation enhancement product is similar to what 
has been produced for many years for land-usage studies with Landsat TM data. The 670 
ran (red) channel in the previously discussed RGB true color enhancement is replaced by 
the 865 nm (SeaWiFS channel 8) channel data. This results in green vegetation 
appearing as bright red and appears usually in sharp contrast to surrounding regions 
devoid of vegetation. An example of this product is shown in Fig. 6. The vegetation 
enhancement may prove useful in identifying potential fetch areas for airborne dust. 
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|0061] Airborne Dust Enhancement 

|0062] An important capability of SeaWiFS, resulting from its multiple narrowband 
visible channels and indicated in the examples above, is discrimination between airborne 
dust and clouds/smoke via true color methods. Fine-grain particles lifted from sparsely 
vegetated and dry desert surfaces by strong frictional winds (typically greater than 10 
m/s, usually by meso- and synoptic-scale circulation features) give rise to significant 
plumes of atmospheric dust. These plumes follow the low to mid-level atmospheric 
circulation and may also become incorporated within synoptic scale baroclinic systems 
(giving rise on occasion to unusually "mud rain" storms that have been reported to 
deposit a coating upon everything in their path). 

[0063] Perhaps more importantly to Navy interests, airborne dust has significant 
detrimental impacts on slant-range visibilities (with visible optical depths that can exceed 
3.0 and accompanied by visibilities less than 1 nautical mile) and also poses a serious 
hazard in terms of its potential to damage turbine engines. This is of particular relevance 
to aircraft carrier operations in regions of the world prone to widespread dust outbreaks. 
Because airborne dust retains the earth tone hues (usually light brown or yellow, 
depending on the source region) of the underlying surface, it is readily identifiable in 
most cases from a true color image, especially for optically thicker dust plumes. For 
regions of less dust (where its effects on slant range visibility may still be appreciable), 
an additional enhancement beyond true color is required. NRL Monterey has developed 
a dust enhancement algorithm for SeaWiFS that attempts to provide this capability. 
Enhancement of dust is predicated on basic spectral properties that characterize its 
appearance at visible wavelengths. As indicated in Fig. 7, yellow light (characteristic of 
many desert dust plumes) results from the removal of blue light (via preferential 
absorption by the dust). This spectral contrast is even more dramatic between the blue 
and shortwave infrared channels. This preferential absorption does not occur to the same 
magnitude in liquid/ice clouds (explaining why clouds appear as shades of gray to white, 
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indicating nearly equal reflected intensities of red/green/blue light, in a true color 
composite). By defining a dust enhancement parameter "A": 




where a>. are reflectances at wavelength X (nm), any dust component of the image will be 
enhanced. Cloud pixel values, which are relatively spectrally flat, will be diminished in 
brightness owing to the small numerator difference. The logarithm of the A parameter is 
scaled between -0.45 and 0.20 and loaded into the red channel of the RGB composite. 

[0064] Owing to the significantly reduced red component for cloud pixels in the dust 
enhancement (small values of the A parameter) image, clouds typically appear as shades 
of light blue or cyan. Retaining the original gray/white cloud hues in the dust 
enhancement (for cosmetic purposes only) requires the implementation of a simple cloud 
mask — pixels flagged as cloudy can be re-assigned to their original red channel values. 
Clouds in the imagery are identified first by computing the mean scaled magnitude and 
standard deviation of the 412. 555, and 670 nm channels. Mean magnitudes exceeding 
50% and standard deviation less than 2.5% are used as thresholds for positive cloud 
identification. The simple premise behind choosing these thresholds is that clouds are in 
general bright (high brightness magnitude) and relatively spectrally flat over the visible 
(low standard deviation). Pixels flagged as "cloudy" are reset to the red channel value in 
the A parameter. As the dust enhancement enhances land surfaces also, a land mask is 
applied to this product. 

[0065] Fig. 8 demonstrates the dust enhancement product covering the Persian Gulf 
region. The top panel pair reveals that while the true color imagery is useful in 
identifying the significant dust plumes, the dust enhancement reveals a far more 
extensive coverage of airborne dust circulating within the Gulf. The center and lower 
panel pairs further illustrate this enhancement capability. While not necessarily posing 
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the same hazard to aircraft engines as the optically thick (bright) plumes, these areas of 
enhanced aerosol have detrimental impacts on slant-range visibility. Although more 
difficult to detect than over the oceans, the presence of significant dust over bright and/or 
sparsely vegetated land surfaces can in some cases be inferred by noting contrast 
reduction in the background terrain. The lack of an infrared channel on the SeaWiFS 
instrument precludes the technique of detecting elevated (cooler) dust above a warm 
surface background via emission measurements. The utility of EOS-MODIS, which 
provides 36 channels (including narrowband red, green blue, NIR and LWIR bands), will 
be explored in this capacity. 

[0066] Other Utilities 

[0067] The global coverage afforded by the polar orbiting platform provides SeaWiFS 
with opportunities to survey many atmospheric and terrestrial phenomena. Its utility as 
an additional sensor for tropical storm monitoring (e.g., Typhoon Pabuk shown in Fig. 9) 
may increase in relevance in light of the current mechanical problems aboard the 
Geostationary Meteorological Satellite (currently providing coverage of the Tropical 
Western Pacific and the primary satellite observation platform supporting the Joint 
Typhoon Warning Center). Unexpected events provide unique opportunities to examine 
additional capabilities of the SeaWiFS channels. The eruption of Mt. Etna, shown in Fig. 
10, demonstrates an additional utility of the dust enhancement product — revealing 
features of the ash plume difficult or impossible to detect in panchromatic or even true 
color renditions. Depending on their composition and concentration (not all volcanic 
plumes are comprised of materials that will be enhanced by the currently detailed dust 
enhancement), such plumes may be tracked around the globe, providing invaluable 
guidance to the commercial airline industry for avoiding these potentially serious flight 
hazards. The detection capabilities of SeaWiFS offer a valuable asset to Volcanic Ash 
Advisory Centers distributed world-wide whose charter is to monitor and track these 
plumes globally, and also the Geneva-based World Meteorological Organization, which 
coordinates global scientific activity among an international panel of participants. 
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[0068] New Processing Formats 

[0069] The SeaWiFS processing as discussed in this report is well suited for the creation 
and analysis of high quality SeaWiFS imagery products in an automated, near real-time 
framework. 

[0070] For the reasons cited above, efforts to reproduce this processing in a fashion 
entirely independent of third party software have resulted in an analogous package that 
assumes receipt of SeaWiFS data in TeraScan Data Format and proceeds to customize the 
data in the FORTRAN-90 programming language. True color and dust enhancement 
image rendering is then completed using the 'Ygbimage" TeraScan intrinsic function. 
This FORTRAN-based version is actually more efficient and hence runs faster than its 
IDL counterpart. 

[0071] Blended Products 

[0072] SeaWiFS is only one of several telemetries available for processing at NRL 
Monterey. It is a straightforward procedure to incorporate data from these other sensors 
within the same image to produce a multi-sensor product. For example, the masked-out 
land areas in the dust enhancement product may be filled in by temporally matched 
geostationary data (e.g.. Meteosat 5 over Bahrain, MeteoSat 7 over Rota, and GMS over 
Yokosuka). Infrared information can be used to enhance regions of elevated dust over 
the land in certain situations. As NRL also has at its disposal analysis and forecast data 
from its mesoscale (Coupled Ocean-Atmosphere Mesoscale Prediction System)™ and 
global scale (Navy Operation Global Atmospheric Prediction System) models, it may 
also prove useful to blend these data into the products. Using again the example of the 
dust enhancement product, inclusion of the COAMPS™ 925 mb wind field in the form of 
vectors overlaid upon the imagery may serve not only as a qualitative validation of the 
model analysis but also provide the METOC officer in the field with a means to 
predicting the migratory behavior of significant dust plumes in the short term. Overlay of 
short-term forecasts of the frictional wind over land (winds near the surface that are 
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responsible for lifting the dust into the atmosphere) may also be of use. Of course, any 
such multi-sensor and model-fusion products will require a training module identifying 
the various components, strengths, and weaknesses. 

[0073] These data are of immediate relevance to both the initialization and validation of 
the Navy Aerosol Analysis and Prediction System (NAAPS; a global, multi-component 
analysis and modeling capability to produce aerosol products from satellite data and 
surface-based measurements). The goal of NAAPS is to improve the 
forecasting/nowcasting of visibility restrictions caused by aerosols, including dust, 
smoke, and volcanic ash. 

1 0074] Extension to MODIS (moderate resolution imaging spectroradiometer) Processing 

[0075] The MODIS instrument offers a similar true color imaging capability together 
with improved spatial resolution (up to a factor of 4 over SeaWiFS HRPT). Fig. 1 1 
provides an example of the MODIS dust enhancement over the Northern Arabian Sea. 
Plumes of dust over ocean and dry lakebeds on the desert interior (often serving as 
sources for dust) are enhanced as shades of pink. Of specific relevance to the dust 
enhancement product is the availability of thermal infrared channels that enable the 
discrimination of elevated (cooler) dust plumes over warm backgrounds. On-board 
calibration of MODIS radiances supports accurate physical retrievals of aerosol optical 
properties. A technique for producing 250m true color imagery using 500m blue/green and 250m 
red channels has also been developed. This high spatial resolution capability is enhanced by 
temporal coverage when considering the Terra (MODIS: 1030 descending node). SeaStar 
(SeaWiFS: 1200 descending node), and Aqua (MODIS: 1330 ascending node; platform to be 
launched in April 2002) platforms in concert. This quasi-geostationary 3-hour loop will provide 
valuable information on cloud and/or dust motion. 

[0076] Conclusion 
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[00771 A new processing capability for generation of true color, vegetation enhancement, 
and dust enhancement imagery from the 8-channel SeaWiFS instrument has been 
developed by the NRL Monterey Satellite Meteorology Applications Section team. 

[0078] Two main points become clear from the above commentary: I) SeaWiFS and 
other sensors offering true color imaging capability are of significant value and in high 
demand for Navy operations, and 2) in spite of the fact that such data has been available 
for some time now from the satellite platform (e.g., LandSat-TM) this capability has not 
been transitioned effectively to the end-users who require it (i.e., those deployed in the 
field). In turn, the twofold intent of this report has been to demonstrate the current in- 
house processing capabilities for the SeaWiFS sensor and articulate this operational 
demand. A logical extension to these statements is that there exists a pressing need for 
X-band receiving stations in perennially sensitive regions of the world for the purpose of 
capturing MODIS and other high-data-rate (both multi- and hyper-spectral) satellite 
telemetries. The potential savings in terms of both loss of life and equipment (e.g.. Cantu 
[2001] reports a 51 million dollars annual loss in Navy equipment attributed to visibility- 
related meteorological phenomena such as aerosols, dust, and atmospheric trapping 
layers) owing to improved METOC guidance gained from these state-of-the-art sensors 
far outweighs the one-time cost of the receiving stations themselves. 



